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I.  INTRODUCTION 


A.  MOTIVATION 

The  start  of  the  fourteenth  century  marked  the  end  of  hand-to-hand  combat  and 
ushered  in  a  new  martial  law,  with  the  introduction  of  the  cannon  proper  on  the  battlefield 
[1].  By  the  seventeenth  century,  artillery  was  employed  with  great  effect  against  enemy 
troops  in  the  open  [1], 

Since  then,  man  has  been  trying  to  improve  the  precision  and  range  of  artillery. 
Modern  warfare,  characterized  by  urban  close  fire  support  and  larger  areas  of  operation, 
has  further  escalated  the  demands  for  precision  and  range  [2],  In  the  current  context,  a 
precision  guided  missile  presents  itself  as  an  effective  but  costly  option  [3].  A  more 
economical  option  is  to  explore  ways  to  improve  the  existing  projectile  stockpile’s 
precision  and  range  by  redesign. 

B.  TRADITIONAL  PROJECTILE  DESIGN 

A  conventional  155mm  projectile  consists  of  a  hollow  steel  shell  filled  with  high 
explosive  [4].  The  projectile’s  body  generally  begins  with  a  streamlined  nose  and 
terminates  using  a  boat  tail  base  for  aerodynamic  efficiency.  Prior  to  firing,  a  fuze  is 
placed  into  the  fuze  adapter  located  at  the  projectile  nose. 

Conventional  155mm  projectile  design  was  popularized  during  the  World  War  II 
era.  Since  then,  modifications  to  the  projectile  design  have  been  minimal.  The  projectile 
is  “dumb”  and  employed  in  an  open  loop  manner.  Once  in  flight,  the  dumb  round  does 
not  have  any  means  to  guide  itself  to  the  target.  Therefore  in  a  typical  engagement 
scenario,  several  iterative  adjustments  are  made  before  the  projectile  is  able  to  land  near 
the  intended  target.  This  inevitably  drives  up  the  mission  cost,  lengthens  the  mission 
duration  and  increases  the  possibility  of  enemy  counter-fire.  This  is  even  more 
unacceptable  in  urbanized  combat  zones,  where  collateral  damage  is  less  tolerable. 
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C.  DESIGN  CONSIDERATIONS 

The  flight  of  a  projectile  is  complex  in  nature.  The  projectile  is  subjected  to  a 
wide  range  of  Mach  numbers  during  its  flight.  This  complexity  is  further  compounded  by 
the  fact  that  the  projectile  is  subjected  to  varying  atmospheric  conditions  as  it  cuts 
through  various  altitudes. 

Therefore,  an  understanding  of  the  aerodynamic  characteristics  of  a  projectile  is 
critical  in  the  design  of  new  projectile  configurations  to  achieve  improved  perfonnance 
[5],  From  studies,  the  critical  aerodynamic  behavior  occurs  in  the  transonic  speed  range 
from  Mach  0.9  to  1.2  [5].  An  in-depth  understanding  of  the  aerodynamic  behavior  in  this 
region  is  crucial  in  the  development  of  a  projectile  [5],  Aerodynamic  properties  of  a 
projectile  have  been  found  to  vary  drastically  within  the  transonic  speed  region  [6].  With 
current  advancements  in  simulation  technologies  and  perhaps  a  more  in  depth 
understanding  of  the  subject,  it  may  be  possible  to  make  a  minor  modification  to  a 
traditional  projectile  to  improve  its  performance,  using  existing  artillery  howitzers. 

D.  PREVIOUS  WORK 

This  thesis  continued  the  aerodynamic  studies  of  the  M549  projectile  started  by 
[5],  [6]  and  [7].  Past  studies  focused  on  finding  the  aerodynamic  properties  of  the 
projectile  using  dated  modeling  simulation  and  modeling  tools.  Current  advancement  in 
simulation  and  modeling  justified  a  reinvestigation. 

Recent  work  of  [7]  focused  solely  on  the  drag  prediction  for  the  M549  projectile. 
A  more  in-depth  understanding  of  other  important  aerodynamic  coefficients  like  pitching 
moment  and  normal  force  is  required.  In  addition,  the  simulation  and  modeling  tools  used 
by  [5]  and  [6]  have  been  superseded  by  newer  versions.  These  observations  justify  further 
investigation. 
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E.  SCOPE  AND  ORGANIZATION 

1.  Scope 

The  scope  of  this  thesis  aims  to  get  a  better  understanding  of  an  M549  projectile’s 
stability  derivatives.  The  modeling  program  SolidWorks  is  used  to  model  the  projectile. 
The  simulation  program  ANSYS  CFX  is  employed  to  detennine  the  stability  derivatives 
of  the  projectile  at  various  conditions.  AP09  and  live  firing  data  are  used  as  benchmarks 
for  the  simulation  results. 

Comparison  studies  are  conducted  between  two  angle  of  attack  (AO A)  generation 
methodologies  in  ANSYS  CFX.  Finally,  an  investigation  is  conducted  to  explore  the 
effect  of  nose  modification  on  the  M549  155mm  projectile’s  lift  coefficient.  The 
redesigned  projectile’s  aerodynamic  data  is  compared  against  the  baseline  projectile’s 
data. 


2.  Organization 

Chapter  II:  The  separation  of  the  aerodynamic  forces  into  lift  and  drag  is 
presented.  This  is  followed  by  an  explanation  of  the  aerodynamic  moment.  Equations 
relevant  to  both  aerodynamic  force  and  moment  are  described. 

Chapter  III:  The  use  of  the  simulation  program,  ANSYS  CFX  in  the  aerodynamic 
prediction  of  the  projectile  is  described.  Advice  that  is  helpful  in  resolving  some  of  the 
common  simulation  issues  is  documented.  In  addition,  optimization  and  tuning 
techniques  to  improve  the  accuracy  of  the  simulation  is  discussed.  Finally,  an 
introduction  of  two  types  of  AOA  generation  techniques  is  presented. 

Chapter  IV:  The  use  of  the  simulation  program  AP09  in  the  aerodynamic 
prediction  of  the  projectile  is  described. 

Chapter  V:  The  results  and  analysis  for  both  the  standard  and  modified  projectile 
configuration  are  discussed.  Results  from  the  two  different  methods  of  AOA  generation 
techniques  are  evaluated.  The  accuracy  of  the  simulation  results  is  validated  by 
comparing  the  results  with  published  experimental  results. 
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Chapter  VI:  The  conclusions  drawn  from  this  research  and  recommendations  for 
follow-on  research  are  discussed. 
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II.  DEVELOPMENT  OF  AERODYNAMICS 


Aerodynamics,  a  branch  of  theoretical  physics,  has  progressed  greatly  since 
Leonardo  da  Vinci  speculated  on  the  possibility  of  human  flight  [8].  Aerodynamics  is 
concerned  with  the  interaction  between  a  moving  object  and  the  surrounding  air. 
Aerodynamics  can  be  further  divided  into  two  main  areas,  external  and  internal 
aerodynamics.  Internal  aerodynamics  focuses  on  the  study  of  flow  though  passages 
within  an  object.  In  this  thesis,  the  concern  is  the  external  aerodynamic  flow  around  a 
projectile. 

Nature  can  only  transmit  a  force  to  a  moving  body  through  the  pressure  and  shear 
stress  distributions  on  the  body  surface.  Regardless  of  the  complexity  of  a  body  shape, 
aerodynamic  forces  and  moments  on  the  body  are  due  to  pressure  and  shear  stress 
distributions  over  the  body  surface.  [9] 

A.  AERODYNAMIC  FORCE 

An  aerodynamic  force  is  the  result  of  the  integrated  pressure  and  shear  stress 
distribution  over  the  entire  body  surface.  The  aerodynamic  force  can  be  divided  into  lift 
and  drag  as  shown  in  Figure  1  [9]. 


N 


Figure  1.  Aerodynamic  Forces.  From  [9] 
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1. 


Lift 


Lift  is  defined  as  the  force  component  which  acts  on  a  body  in  a  direction 
perpendicular  to  the  free-stream  velocity  vector.  Regardless  of  the  variation  in  AOA,  lift 
always  maintains  this  orientation  to  the  velocity  vector  [9],  The  lift  coefficient,  Cl  is  a 
dimensionless  coefficient  that  relates  lift  generated  by  a  body  highlighted  in  Equation  2.1. 


CL  = 


0-5pV/Aref 


(2.1) 


2.  Drag 

In  a  real  viscous  flow,  total  drag  consists  of  friction  and  pressure  drag  force  [10]. 
Pressure  drag  is  created  by  normal  force  perpendicular  to  the  boundary  surface  [8]. 
Friction  drag  is  created  by  tangential  forces  parallel  to  the  boundary  surface  [8].  In  a 
viscous  flow,  particles  come  to  a  stop  at  the  body  surface.  The  sum  of  this  effect 
generates  the  skin  friction  drag.  At  certain  conditions,  the  flow  is  separated  from  the  body 
surface  creating  a  pressure  difference  which  causes  the  pressure  drag  [10].  It  is  important 
to  highlight  that  even  when  the  flow  is  not  detached  from  the  body  surface,  shear  forces 
can  still  alter  the  pressure  difference  and  create  pressure  drag  [10].  Drag  is  defined  as  the 
force  component  which  acts  on  a  body  in  a  direction  parallel  to  the  free-stream  velocity 
vector  [9].  Similar  to  lift  force,  regardless  of  the  variation  in  AOA,  drag  always 
maintains  this  orientation  to  the  velocity  vector. 


CD 


D 

0-5  pV>ref 


(2.2) 


3.  Normal  and  Axial  Force 

Normal  force  can  be  defined  as  the  force  component  which  acts  on  a  body  in  a 
direction  perpendicular  to  the  body  axis.  Axial  force  is  defined  as  the  force  component 
which  acts  on  a  body  in  a  direction  parallel  to  the  axis  line.  Regardless  of  the  variation  in 
AOA,  both  nonnal  and  axial  force  always  maintains  this  orientation  to  the  axis  line.  The 
following  equations  relate  the  lift  and  drag  force  to  the  normal  and  axial  force. 


Lift  =  Fn  cos(a)-FA  sin(a) 
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(2.3) 


Drag  =  FN  sin(a)  +  FA  cos(or) 


(2.4) 


B.  SLENDER  BODY  THEORY 

For  slender  bodies  of  revolution  flying  at  small  AOA,  the  lift  can  be  estimated  by 
considering  only  the  flow  in  a  plane  perpendicular  to  the  flight  direction,  the  cross  flow 
plane  [10].  This  reduces  the  problem  to  finding  the  solution  for  incompressible  flow  over 
a  cylinder.  Using  slender  body  theory  in  Figure  2,  the  lift  coefficient  is  found  to  be 
linearly  proportional  to  the  AOA  [10]. 


View  of  d(x) 


Figure  2.  Slender  Body. 

From  [10],  the  derivation  of  lift  coefficient  is  given  by: 

Local  lift  is  defined  as  L'  =  ^L[x) 

Sx 

=  Jq  pR(x) cos  Odd 
=  ^(P-  PMX)  cos  OdO 


Since,  Cp=f^- 

~PyJ 


)R(x)c os  6d6 
(Cp  )R(x)c os  Odd 


(2.5) 


(2.6) 

(2.7) 


(2.8) 
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From  the  incompressible  flow  over  a  cylinder,  one  finds  that, 


C(r  =  R(x))  = 


2  cos  6  8  2 


VnR(x)  8x 


(R  (x)wa(x)) 


L'(x>=^r 


2  Jo 


2  cos  0  8  ,^2 


VxR(x)  8x 


(R  (x)wa(x)) 


R(x)  cos  9d6 


r2  k 


cos2  8d8-^—R2{x)w  lx)) 

8x 


(2.9) 

(2.10) 


=  pV*x-?-R2(x)wa(x) 

8x 


where  for  flight  at  a  steady  AO  A,  wa  =  Vjx  , 


r  Tail  8L  _ 

L  =  — 8x 

J  Nose 


c Tail  8  r, 

=  pVjtl  —R2(x)Vxadx 

J  Nose  sjx 


For  a  body  with  a  pointed  nose  and  a  finite  base  radius  RB, 
L  =  pV2nR2a 


C-T 


pV2xRB2 


pVa2xRB2a 


1 


pVjnRB2 


CL  =2  a 


(2.11) 


(2.12) 


(2.13) 

(2.14) 


(2.15) 


Equation  2.15  provides  a  quick  estimate  of  the  lift  at  small  AOA  [9,  10].  However, 
it  predicts  no  effect  of  Mach  number  and  Reynolds  number  on  the  body  geometry. 
Therefore,  more  accurate  theories  are  needed  to  account  for  these  effects. 
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c. 


AERODYNAMIC  MOMENT 


A  force  has  both  direction  and  magnitude.  The  product  of  force  and  perpendicular 
distance  from  the  center  of  gravity  is  known  as  moment  or  torque.  The  aerodynamic 
moment  which  act  to  change  the  AOA  is  called  a  pitching  moment  [9] .  Positive  pitching 
moment  occurs  when  a  body  is  pitched  in  a  nose  up  direction.  The  pitching  moment 
coefficient  then  is: 


C,.  — 


M 


0-5pVx  ArefLref 


(2.16) 


D.  COMPRESSIBILITY  EFFECT 

The  disturbances  emitted  by  a  disturbance  source  travel  at  the  speed  of  sound,  as 
illustrated  in  Figure  3  [11,  12].  The  generated  wave  patterns  are  strongly  dependent  on 
the  ratio  of  the  speed  of  the  disturbance  source  to  the  speed  of  sound.  With  that,  it  is 
important  to  introduce  the  definition  of  Mach  number.  It  is  named  after  Ernst  Mach  and  is 
defined  as  the  ratio  between  the  velocity  of  a  body  moving  through  the  air  and  the 
velocity  of  sound  in  the  air  [8]. 

M  =  —  (2.17) 

a 

As  depicted  in  Figure  3,  if  M<1  or  V<a,  the  disturbances  can  propagate  ahead  of 
the  moving  body  whereas  at  M=1  the  body  moves  as  fast  as  the  disturbance  fronts.  If 
M>1  or  V>a,  the  disturbances  can  only  propagate  downstream  of  the  body  within  the  so- 
called  Mach  cone.  It  is  obvious  that  fundamentally  different  flow  phenomena  occur 
depending  on  the  value  of  the  Mach  number,  requiring  more  accurate  analyses  than 
slender  body  theory.  Such  analyses  have  been  developed  in  recent  years  which  can  not 
only  predict  the  Mach  number  effects  but  also  viscous  flow  effects.  They  are  based  on  the 
numerical  solution  of  the  Navier-Stokes  equations  for  compressible  viscous  flows  [10]. 
To  this  end,  a  system  of  six  equations  consisting  of  the  continuity  equations,  the  three 
momentum  equations,  the  energy  equations  and  the  equation  of  state  for  air  needs  to  be 
solved  [10], 
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Continuity  Equation: 


dp  |  d(pu)  |  d(pv)  |  d(pw)  =  Q 
dt  dx  dy  & 


Momentum  Equation: 

du  du  du  du  dP 

p - b  pu - b  pv - b  pw — = - 

dt  dx  dy  dz  dx 

d  f  du  ( 2\  f du  dv  dw3] 

dx  y  dx  ^  3 J  ^  dx  dy  dz  J  J 

d  f  f  du  dv3]  d  \  ( dw  dw3] 

dy  {  l^dy  dx  J  J  dz  [  ^dx  dz  )J 

dv  dv  dv  dv  dP 

p - b  pu - b  pv - b  pw — = - 

dt  dx  dy  dz  dy 

d  f  dv  ( 2^)  f du  dv  dw3] 

dy  \  dy  Uj  ^dx  dy  dz  J] 

d  f  f du  dv3]  d  f  f dw  dv3] 

dx  y  ^dy  dx  J  J  dz  y  ^dy  dz  J  J 

dw  dw  dw  dw  dP 

P - b  PU - V  PV - b  PW - = - 

dt  dx  dy  dz  dz 


23  f  du  dv  dw3] 
3  J  ^dx  dy  dz  J  J 


d  \  ( dw  du\\  8  \  f dv  dw3] 
dx  {  ^  dx  dz  J]  dy  y  ^dz  dy  J  J 


(2.18) 


(2.19) 
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Energy  Equation: 


(2.20) 

(2.21) 
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M>1 


Figure  3.  Subsonic,  Sonic  and  Supersonic  Motions.  From  [11] 
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III.  ANSYS  CFX 


One  of  the  commercially  available  codes  for  the  numerical  solution  of  the  system 
of  six  equations  is  the  ANSYS  CFX  code. 

Modern  Computational  Fluid  Dynamics  (CFD)  can  be  regarded  as  equal  partner 
with  pure  theory  and  pure  experiment  in  the  analysis  of  aerodynamic  flows.  As  illustrated 
in  Figure  4,  they  do  not  stand  alone  but  help  each  other  to  understand  aerodynamic  flows 
[9]. 

Pure  Pure 

Experiment  Theory 

\  y 

CFD 

Figure  4.  Aerodynamic  Partnership.  From  [9] 

A.  WORK  FLOW 

ANSYS  CFX  is  a  CFD  solver  program  that  utilizes  the  finite  volume  technique 
[13],  [14].  In  this  technique,  the  area  of  interest  is  subdivided  into  smaller  regions.  Each 
region  is  discretized  and  iteratively  solved.  In  this  manner,  an  approximation  of  the  value 
of  each  variable  at  a  specific  point  within  the  domain  can  be  detennined.  One  can  then 
develop  a  full  understanding  of  the  flow  behavior.  While  the  system  of  Equations  2.18- 
2.21  completely  describes  laminar  flows,  the  analysis  of  turbulent  flow  requires  an 
averaging  process  known  as  Reynolds  averaging,  which  introduces  the  need  for 
turbulence  modeling  [15],  The  work  process  flow  of  ANSYS  CFX  is  depicted  in  Figure 
5. 
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Figure  5.  ANSYS  CFX  Work  Process  Flow.  From  [14] 
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1. 


Model  Generation 


In  this  thesis,  the  body  of  interest  is  a  modified  M549  155mm  projectile  with  a 
power  series  blunt  nose  and  removed  rotating  band  [16].  The  power  series  blunt  nose 
profile  provides  an  accurate  description  of  the  discontinuities  between  fuze  and  projectile. 
The  entire  projectile  has  a  length  of  874.975mm  and  a  diameter  of  155mm.  The  reference 
area  is  defined  as  the  circular  cross  sectional  area  of  the  projectile.  The  center  of  gravity 
(CG)  of  the  projectile  is  542.5mm  from  the  projectile  nose.  A  diagram  of  the  actual  M549 
projectile  is  depicted  in  Figure  6. 


Figure  6.  Geometry  Profile  of  M549  Projectile.  From  [16] 


The  body  of  interest  is  modeled  using  the  modeling  tool  SolidWorks  before  it  is 
imported  to  ANSYS  CFX  Geometry  for  further  processing.  SolidWorks  is  a  widely  used 
modeler  program  that  applies  a  parametric  feature-based  approach  in  models  and 
assemblies  creation  [17],  [18].  The  model’s  center  of  origin  is  centered  at  the  CG  of  the 
projectile.  The  constructed  model  in  SolidWorks  is  depicted  in  Figure  7. 
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Figure  7.  Constructed  Model  of  Standard  M549  Projectile  in  SolidWorks. 
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The  constructed  model  is  saved  and  imported  into  ANSYS  CFX  Geometry.  It  is 
recommended  to  save  the  model  as  “Parasolid  (*.x_t)”  file  type  to  prevent  compatibility 
issues  during  the  transfer  from  SolidWorks  to  ANSYS  CFX. 

A  control  volume  for  the  projectile  is  also  created  using  SolidWorks.  A  control 
volume  is  fixed  in  space  with  air  moving  through  it  [9].  The  fundamental  physical 
principles  and  required  flow  directions  are  applied  to  the  air  crossing  the  control  surface. 
It  is  recommended  to  have  a  large  control  volume  to  effectively  capture  the  shock  wave 
phenomena.  The  control  volume  is  modeled  to  be  approximately  twenty  times  the  size  of 
the  projectile.  Likewise,  the  model  is  saved  as  “Parasolid  (*.x_t)”  file  type  and  imported 
into  ANSYS  CFX  Geometry.  Both  projectile  and  control  volume  are  then  merged  and  cut 
in  the  ANSYS  CFX  Geometry.  This  effectively  creates  a  2D  slice  through  the  combined 
model  to  reduce  computational  time.  However,  this  imposes  a  need  to  multiply  result 
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values  by  a  factor  of  two  to  properly  account  for  full  size  model.  The  “sliced”  projectile 
within  the  control  volume  is  depicted  in  Figure  8. 
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Figure  8.  Control  Volume  with  “Sliced”  Projectile  in  ANSYS  CFX  Geometry. 


2.  Mesh  Generation 

Mesh  generation  of  the  combined  model  was  performed  using  ANSYS  CFX 
Mesh.  A  mesh  is  a  generated  collection  of  edges,  vertices  and  faces  that  define  the  shape 
and  size  of  the  model.  A  more  detailed  mesh  of  the  model  allows  a  more  accurate 
representation  of  the  model  in  simulation.  A  poor  mesh  often  results  in  abrupt  simulation 
termination  and  poor  convergence.  The  average  generated  mesh  consists  of 
approximately  950,000  nodes  and  5,300,000  elements  to  achieve  an  average  Y+  of  1. 

Additional  refinement  settings  like  Inflation  and  Vertex  sizing  were  applied  to 
specific  areas  on  the  projectile  to  improve  the  number  of  mesh  elements  and  nodes  in  that 
particular  area  of  interest.  The  meshed  profile  of  the  combined  model  is  depicted  in 
Figure  9. 
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Figure  9.  Generated  Mesh  Profile  of  Combined  Model  in  ANSYS  CFX  Mesh. 


3.  Pre-Processor  Setup 

A  simulation  profile  is  created  by  specifying  the  required  boundary  conditions, 
flow  physics,  initial  values  and  solver  parameter  in  ANSYS  CFX  Pre.  In  the  Domain 
Setup  under  Basic  Settings  tab,  Material  is  set  to  Air  Ideal  Gas.  Reference  Pressure  is  set 
to  1  atm.  In  the  Fluid  Models  tab,  Heat  Transfer  Option  is  set  to  Total  Energy  for 
compressible  simulations.  Viscous  Work  Term  is  enabled  to  include  viscous  heating 
effects.  High  Speed  Models  is  activated  in  Expert  parameter  with  a  maximum  continuity 
loop  of  value  3,  Turbulence  Option  is  set  to  Shear  Stress  Transport  (SST).  The  SST 
model  is  known  for  robust  simulation  of  aerodynamic  flows.  It  works  especially  well  to 
predict  separated  flow  behavior.  The  SST  model  effectively  integrates  the  accurate 
fonnulation  of  k-co  model  in  the  near  wall  region  with  the  free  stream  independence  of 
the  k-e  model  in  the  far  field  [19].  This  blending  ensures  that  the  model  equations  behave 
appropriately  in  both  the  near  wall  and  far  field  regions  [20].  With  the  SST  turbulence 
model,  the  simulation  can  handle  any  sensible  value  of  Y+  and  automatically  transition 
from  integration  to  the  wall  to  wall  functions.  Y+  is  easily  extended  to  include  curvature 
effects,  transition  and  several  effects  if  required  [14]. 
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In  order  to  specify  the  flow  direction  and  AO  A,  boundary  conditions  were  created 
under  Fluid  setting.  Under  the  Boundary  Details  tab,  Mass  and  Momentum  option  was 
changed  to  Cart.  Vel  Components  to  simulate  flow  in  u,  v  and  w  directions  were  set  as 
illustrated  in  Figure  10.  Static  temperature  was  set  at  288. 15k. 
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Figure  10.  Boundary  Condition  Setting  in  ANSYS  CFX  Pre. 


4.  Solver  Control 

Upon  completion  of  the  pre-processing  setup,  solver  control  settings  are  initiated 
in  ANSYS  CFX  Pre  for  the  simulation  runs.  The  solver  resolves  the  Navier-Stokes 
equations  using  the  created  finite  volume  and  algebraic  multi-grid  approach  [14].  It 
applies  the  law  of  conservation  of  mass,  momentum  and  energy  to  the  created  control 
volume  [14].  Partial  differential  equations  are  integrated  over  the  control  volumes.  These 
equations  are  translated  to  algebraic  equations  by  approximations.  These  equations  are 
then  iteratively  solved.  The  solver  will  end  the  simulation  and  generate  results  when 
residuals  converge  to  a  steady  state.  Convergence  occurs  when  the  residuals  have  been 
reduced  by  at  least  three  orders  of  magnitude  [14].  The  maximum  iteration  is  set  to  100  to 
ensure  convergence.  The  timescale  factor  is  increased  from  1  to  10  to  accelerate 
convergence  to  steady  state  as  fast  as  possible.  Residual  target  is  set  to  le-7  to  ensure 
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solution  converged  well.  Often,  one  has  to  strike  a  balance  between  iteration  counts, 
residual  target  and  timescale  factor  for  solution  optimization. 

5.  Solver  Manager 

The  flow  computation  is  controlled  in  ANSYS  CFX  Solver  Manager.  The  Solver 
Manager  will  attempt  to  resolve  all  the  solution  variables  for  the  problem  specified  in  the 
Solver  Control.  In  Solver  Manager,  users  can  choose  either  “Initial  Conditions”  or 
“Current  Solution  Data”  under  the  Initialization  Option.  “Current  Solution  Data”  option 
is  recommended  if  a  previous  run  of  the  same  profile  was  conducted.  The  solver  will 
make  use  of  the  past  generated  data  to  kick-start  the  new  simulation.  “Double  Precision” 
is  chosen  to  maintain  simulation  accuracy. 

Run  mode  is  dependent  on  the  number  of  processors.  One  can  select  parallel 
processing  if  more  than  one  processor  is  available  [14].  Increasing  the  number  of 
processors  and  Random  Access  Memory  (RAM)  capacity  will  drastically  reduce  the 
required  computational  time.  The  recommended  computer  hardware  is  at  least  8GB  RAM 
and  four  core  processors.  During  the  simulation,  it  is  recommended  to  open  the  Display 
Monitor  to  monitor  the  progress  of  the  runs  and  detect  any  abnormal  report.  A  snapshot 
of  ANSYS  CFX  Solver  Manager  is  depicted  in  Figure  11.  Upon  the  completion  of  the 
simulation,  the  results  files  are  generated. 


^5  AJ :  c7i/»OAO£  0,  ’  Hunt  ffc*  IC^XI  C*X  M.-mjr 
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Figure  1 1 .  “Run  Define”  in  ANSYS  CFX  Solver  Manager. 
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6. 


Post  Processing 


Upon  the  successful  completion  of  the  simulation  run,  output  results  and  plots  can 
be  seen  from  ANSYS  CFX  Post.  ANSYS  CFX  Post  is  an  interactive  post  processing  tool 
that  allows  the  visualization  of  the  output  results.  Scalar  variables  like  pressure,  velocity, 
Mach  numbers  and  Y+  value  can  be  checked.  A  quantitative  numerical  calculator  known 
as  Function  Calculator  is  used  to  calculate  the  various  flow  parameters.  Interactive  view 
of  the  shockwave  phenomena  can  be  called  up.  The  post  processing  tools  in  ANSYS 
CFX  Post  are  depicted  in  Figure  12. 
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Figure  12.  Post  Processing  Tools  in  ANSYS  CFX  Post. 


21 


7. 


Coordinate  Axis  Definition 


Similar  coordinate  axes  illustrated  in  Figure  13,  are  defined  for  both  SolidWorks 
and  ANSYS  CFX.  The  center  of  gravity  is  defined  at  the  center  of  origin.  The  simulated 
free-stream  flow  is  running  from  left  to  right  for  AOA  of  0°.  Drag  force  or  force 
component  in  -X  direction  can  be  extracted  using  “force  x”  in  the  Function  Calculator. 
Lift  force  or  force  component  in  Y  direction  can  be  extracted  using  “force  y”  in  the 
Function  Calculator.  Torque  about  Z  axis  (into  the  diagram)  can  be  extracted  using 
“torque  z”  in  the  Function  Calculator. 
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B.  OPTIMIZATION  AND  TUNING  TECHNIQUES 


The  state  of  a  model  mesh  has  a  direct  effect  on  the  generated  results.  A  high 
resolution  mesh  is  critical  to  the  integrity  of  the  generated  output.  In  order  to  improve  the 
accuracy,  one  should  aim  to  enhance  meshing  within  computational  and  resource 
constraint.  However,  a  general  enhancement  of  mesh  will  vastly  increase  computational 
load.  Therefore  it  is  advisable  to  enhance  meshing  at  selected  areas  of  interest  that  might 
experience  drastic  variations  in  pressure  or  velocity.  A  coarser  mesh  may  be  employed  on 
other  non-critical  areas  to  reduce  computational  time.  Several  mesh  refinement 
techniques  were  tested  and  two  techniques  stood  out  in  terms  of  performance  and  ease  of 
employment. 

1.  Inflation 

Inflation  is  a  useful  mesh  enhancement  technique  for  improving  boundary  layer 
resolution  [14].  Inflation  improves  the  mesh  of  the  boundary  layer  by  further  refining  the 
mesh  normal  to  the  wall.  Inflation  can  be  employed  after  the  initial  meshing  of  the  model 
is  completed.  It  is  especially  useful  in  reducing  the  Y+  value  to  an  acceptable  level  in 
order  to  resolve  the  viscous  sub-layer.  Y+  is  the  dimensionless  distance  from  the  wall  and 
it  is  an  indicator  of  how  well  the  boundary  layer  is  being  modeled.  Y+  should  be 
approximately  1  or  below  to  capture  laminar  and  transitional  boundary  layers  correctly 
[14].  It  is  also  necessary  for  flow  characterized  by  separation  and  reattachment.  From 
Equation  3.1,  Y+  is  dependent  on  friction  velocity  [14].  An  increase  in  Mach  number  will 
lead  to  an  increase  in  friction  velocity.  It  is  recommended  to  lower  Y+  from  onstart  with 
a  good  mesh  of  the  model  with  sufficient  inflation  layers. 

(3.D 

V 

It  is  recommended  to  slowly  increase  the  inflation  level  from  default  number.  An 
abrupt  jump  in  the  inflation  layer  number  will  drastically  increase  computational  time.  A 
close  up  view  of  the  generated  inflation  layers  at  boundary  is  depicted  in  Figure  14. 
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Figure  14. 


Inflation  Layers  at  Boundary. 


2.  Sizing 

ANSYS  CFX  Meshing  allows  two  types  of  size  functions  to  provide  appropriate 
mesh  sizing  for  different  physics  [14].  The  default  size  function  is  designed  to  accurately 
capture  the  geometry  while  minimizing  the  number  of  elements  in  the  model  [14].  User 
invoked  advanced  size  function  is  similar  except  that  a  smooth  growth  rate  is  maintained 
between  the  regions  of  curvature  [14].  Region  of  interest  is  determined  by  user  and  the 
mesh  in  the  region  can  be  further  enhanced  with  detailed  meshing.  An  illustration  of 
vertex  sizing  is  depicted  in  Figure  15  and  enhanced  meshing  in  area  of  interest  is  depicted 
in  Figure  16. 
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Figure  15.  Vertex  Sizing 
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Figure  16.  Enhanced  Mesh  at  Area  of  Interest. 
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c. 


ANGLE  OF  ATTACK  GENERATION  METHODOLOGIES 


Generally,  the  body  of  interest  can  be  made  to  experience  an  AOA  by  two 
methods  in  ANSYS  CFX.  The  first  method  is  to  rotate  the  body  of  interest  about  a 
particular  axis.  The  other  method  is  to  modify  the  inlet  velocity  components. 


1.  Body  of  Interest  Rotation 


Using  this  method,  the  free-stream  flow  direction  is  maintained  at  X  axis  direction. 
The  body  of  interest  is  made  to  rotate  to  the  desired  AOA.  Under  ANSYS  CFX  Geometry, 
user  can  insert  a  plane  and  indicate  an  axis  of  interest  for  rotation.  The  degree  of  rotation 
is  the  desired  degree  of  AOA.  This  is  relatively  easy  to  implement.  However,  each 
modification  of  the  rotation  will  require  a  new  mesh.  An  illustration  of  body  rotation  is 
depicted  in  Figure  17. 
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Figure  17.  Body  of  Interest  at  10°  Rotation. 
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2. 


Inlet  Velocity  Modification 


Using  this  method,  the  free-stream  flow  is  deliberately  offset  from  the  X  axis 
direction  to  simulate  the  AOA.  User  will  not  be  required  to  make  any  changes  to  the 
existing  generated  geometry  and  mesh  of  the  body  of  interest.  The  only  required  change 
is  made  in  ANSYS  CFX  Setup.  Under  all  inlet  and  initialization  details,  the  user  is 
required  to  insert  both  u  and  v  velocity  components.  For  a  certain  AOA  profile,  one  can 
make  use  of  the  following  equations  to  detennine  the  required  u  and  v  velocity 
components. 


U  =  Va  cos  (a) 

(3.2) 

V  =  Vn  sin  (a) 

(3.3) 

The  corresponding  creation  of  inlet  is  required  to  match  with  an  outlet.  Therefore 
with  two  inlets,  one  will  expect  two  outlets  to  be  created.  The  u  and  v  inlets  and  outlets 
are  illustrated  in  Figure  18. 
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Figure  18.  U  and  V  Inlets  and  Outlets. 
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IV.  AERO-PREDICTION  (AP)  09  CODE 


The  Aero-Prediction  (AP)  code  is  used  for  the  aerodynamic  prediction  of  mortars, 
low  drag  bombs,  projectile  and  missile  bodies  [21].  The  AP  code  has  undergone  several 
iterations  of  improvement  and  the  thesis  will  make  use  of  the  latest  version  known  as 
AP09  [21].  It  employs  semi-empirical  techniques  based  on  both  experimental  results  and 
theoretical  methods  to  predict  the  aerodynamic  coefficients  [21].  The  greatest  advantage 
offered  by  the  AP09  code  is  rapid  result  generation.  After  the  geometry  and  essential 
simulation  conditions  are  entered  into  the  program,  the  result  is  generated  instantly. 

A.  LOGIC  FLOW 

The  AP09  code  consists  of  the  pre-processor,  the  post-processor,  the  aerodynamic 
module,  trim  aerodynamics  module,  a  ballistic  trajectory  module  and  a  three  degrees  of 
freedom  trim  perfonnance  module  [22].  The  geometry  and  simulation  conditions  are 
keyed  into  the  pre-processor.  The  pre-processor  consists  of  geometry,  free-stream,  option 
and  trajectory  inputs  [22].  After  the  results  generation  is  completed,  the  user  can  access 
the  results  in  the  post-processor.  The  post-processor  consists  of  output  data  which  can  be 
in  the  form  of  tables  of  data  or  plots  of  aerodynamic  data  or  trajectory  information  [22]. 
The  aerodynamics  module  contains  all  aerodynamic  computations  and  is  used  in  a  stand¬ 
alone  mode  or  to  provide  inputs  to  the  trajectory  modules  [22].  User  is  required  to 
generate  a  feasible  geometry  prior  to  result  generation.  The  geometry  inputs  are  required 
for  all  aerodynamic  computations  in  the  code.  AP09  followed  its  predecessor,  AP05’s 
logic  flow  [21].  The  AP  code’s  logic  flow  is  depicted  in  Figure  19. 
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Figure  19.  AP  Code  Logic  Flow.  From  [22] 


B.  PRE  PROCESSOR 

Prior  to  the  creation  of  a  M549  projectile  geometry  model,  the  global 
configuration  is  determined  in  AP09  Pre  Processor.  Under  Configuration  tab  on  the 
AP09’s  interface,  “Body- Alone”  is  selected  for  M549  geometry  modeling.  In  the  Input 
tab,  “Geometry  (Millimeters)”  fonnat  was  selected  as  default  unit  for  the  subsequent  data 
inputs. 
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1.  Geometry  Creation 

The  projectile  dimension  highlighted  in  Figure  3  is  entered  into  the  Body- Alone 
Geometry  Tab.  The  interface  of  Body-Alone  Geometry  Tab  is  illustrated  in  Figure  20. 
With  the  dimension,  AP09  creates  a  geometry  model  of  the  projectile  which  can  be 
viewed  under  Generation  Tab. 
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Figure  20.  Body-Alone  Geometry  Tab. 


Upon  completion  of  Body-Alone  Geometry  Tab’s  entries,  the  Nose  Geometry’s 
entries  are  required  for  the  creation  of  the  projectile’s  nose  profile.  Circular  shape  and 
power  series  blunt  was  selected  for  the  nose  profile.  The  circular  radius  of  the  nose  end  is 
79mm,  half  of  the  projectile’s  diameter.  The  interface  of  Nose  Geometry  is  illustrated  in 
Figure  21. 
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Figure  2 1 .  Nose  Geometry. 


Within  the  nose  profile,  the  power  series  exponent,  truncation  tip  radius  and 
length  of  nose  is  specified.  This  is  illustrated  in  Figure  22. 


Figure  22.  Nose  Profile. 
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Upon  the  completion  of  Nose  Geometry’s  entries,  circular  cross  section  and 
standard  profile  was  selected  for  the  model  as  illustrated  in  Figure  23. 


Figure  23.  Afterbody  Geometry. 


Under  the  Standard  Tab  in  the  Afterbody  Geometry  Tab,  the  longitudinal 
Afterbody  coordinate  from  nose  tip  or  entire  length  of  the  projectile  nose  and  Afterbody 
only  was  indicated  as  illustrated  in  Figure  24. 


Data  is  to  be  entered  in  millimeters  (Format  xxxxxxxx.xxxxxxx) 


Figure  24.  Afterbody  Standard  Tab. 
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Finally,  the  longitudinal  boattail/flare  coordinate  from  nose  tip  or  entire  length  of 
projectile  is  indicated.  The  boattail/flare  half  width  or  radius  is  required  too.  This  is 
illustrated  in  Figure  25. 


Figure  25.  Boattail/Flare  Tab. 

The  essential  dimensions  of  the  projectile  highlighted  in  Figure  3  were  keyed  into 
the  body  alone  geometry  tab.  Using  the  dimensions  given,  AP09  will  generate  a  geometry 
model  which  can  be  viewed  under  Generation  Tab. 

2.  Free  Stream  Condition  and  Option  Setup 

Upon  verification  of  the  geometry  model,  the  free  stream  condition  for  type  of 
simulation  is  specified.  In  this  thesis,  Mach  Sweep  is  initialized  for  a  sweep  from  Mach 
0.7  to  Mach  1.4  at  various  AOA  at  an  altitude  of  20000  feet.  The  simulation  option  is  the 
last  step  before  simulation  generation.  These  are  illustrated  in  Figure  26. 
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Figure  26.  Simulation  Option. 


C.  POST  PROCESSOR 


Upon  the  completion  of  the  result  generation,  the  relevant  result  in  plot  and  table 


format  can  be  called  up.  The  results  can  be  extracted  out  in  text  file  for  plotting  in  Excel. 


35 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


36 


V.  RESULTS  AND  ANALYSIS 


The  generated  results  and  analysis  are  categorized  into  three  main  segments,  Part 
I  for  Standard  Nose  Configuration,  Part  II  for  Comparison  Studies  of  the  AOA 
Generation  Methodologies  and  Part  III  for  Modified  Nose  Configuration.  Experimental 
data  from  [5,  6]  are  used  as  the  benchmark  reference  data  for  the  ANSYS  CFX  and  AP09 
results. 

A.  PART  I-STANDARD  NOSE  CONFIGURATION 

Under  the  Standard  Nose  Configuration,  ANSYS  CFX  and  AP09  simulations 
were  performed  on  a  standard  M549  projectile  model.  Total  drag  coefficient,  normal 
force  coefficient  and  pitching  moment  coefficient  were  evaluated  for  various  AOA  and 
Mach  numbers.  These  simulations  were  performed  using  IVM  method. 

1.  Total  Drag  Coefficient  versus  Mach  Number 

The  generated  total  drag  coefficient  versus  Mach  Numbers  at  AOA  of  0°  results 
are  illustrated  in  Figure  26.  The  variation  of  the  total  drag  coefficient  was  detennined 
using  six  data  points  from  Mach  0.7  to  Mach  1.2  with  an  interval  of  0.1.  The  generated 
Y+  ranged  from  0.8  to  1.6. 
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Figure  27.  Total  Drag  Coefficient  versus  Mach  at  AOA  of  0°. 


From  Figure  27,  ANSYS  CFX  and  AP09  and  reference  data  from  [5]  exhibited 
similar  total  drag  coefficient  trends.  ANSYS  CFX  and  AP09  correctly  predicted  the 
increase  in  total  drag  coefficient  in  transonic  region.  This  increase  can  be  explained  by 
the  change  in  pressure  distribution  caused  by  the  changes  in  flow  speed.  As  the  Mach 
number  reached  transonic  and  supersonic  values,  shocks  start  to  form.  Thus  drag  starts  to 
increase  at  high  speed  flow.  This  form  of  drag  is  known  as  wave  drag. 

Comparing  prediction  accuracy,  ANSYS  CFX  fared  better  in  the  subsonic  and 
transonic  region.  AP09  was  able  to  give  a  better  prediction  in  the  supersonic  region  as 
compared  to  ANSYS  CFX.  This  could  be  due  to  increasing  Y+  value  illustrated  in  Figure 
28. 


38 


Based  on  Equation  3.1,  Y+  is  dependent  on  friction  velocity.  Thus,  increase  in 
Mach  number  will  lead  to  an  increase  in  Y+,  contributing  to  simulation  error.  This  may 
explain  the  large  deviation  of  the  ANSYS  CFX  prediction  from  the  experiment  in  low 
supersonic  flow. 

In  order  to  investigate  whether  these  observations  are  still  valid  at  finite  AOA,  the 
simulations  were  repeated  for  AOA  of  2°  and  4°  as  illustrated  in  Figures  29  and  30.  It 
was  observed  that  the  ANSYS  CFX  and  AP09  prediction  agree  well  at  subsonic  Mach 
numbers,  but  deviate  in  supersonic  flow. 
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MACH 


Figure  29.  Total  Drag  Coefficient  versus  Mach  at  AO  A  of  2°. 
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Figure  30.  Total  Drag  Coefficient  versus  Mach  at  AO  A  of  4°. 
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2. 


Normal  Force  Coefficient  Slope  versus  Mach  Number 


The  normal  force  coefficient  slope  versus  Mach  Number  results  are  illustrated  in 
Figure  31.  The  normal  force  coefficient  slope  was  calculated  by  taking  the  normal  force 
difference  between  AOA  of  2°  and  0°.  The  variation  of  the  normal  force  coefficient  slope 
with  Mach  number  was  determined  using  eight  data  points  from  Mach  0.7  to  Mach  1.4 
with  an  interval  of  0. 1 .  The  simulation  was  conducted  using  ANSYS  CFX  and  AP09.  The 
generated  Y+  ranged  from  1  to  1.8. 
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Figure  3 1 .  Normal  Force  Coefficient  Slope  versus  Mach  Number. 


From  Figure  31,  it  is  apparent  that  the  ANSYS  CFX  and  AP09  predictions  differ 
quite  substantially  from  the  experiment  [5].  However,  ANSYS  CFX  correctly  predicted 
the  decrease  in  normal  force  coefficient  in  the  transonic  region  found  in  the  experiment. 
The  ANSYS  CFX  deviation  from  both  AP09  and  the  experiment  at  Mach=1.4  may  be 
due  to  the  same  reason  explained  before. 
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3.  Pitching  Moment  Coefficient  Slope  versus  Mach  Number 

The  pitching  moment  coefficient  slope  versus  Mach  Number  results  are  illustrated 
in  Figure  32.  The  pitching  moment  coefficient  slope  was  calculated  by  taking  the  torque 
(Z  axis)  difference  between  AOA  of  2°  and  0°.  The  variation  of  the  pitching  moment 
coefficient  slope  from  subsonic  to  supersonic  region  was  determined  using  eight  data 
points  from  Mach  0.7  to  Mach  1.4  with  an  interval  of  0.1.  The  generated  Y+  ranged  from 
1  to  1.8. 


Figure  32.  Pitching  Moment  Coefficient  Slope  versus  Mach  Number. 


From  Figure  32,  it  is  seen  that  ANSYS  CFX  correctly  predicted  the  increase  in 
pitching  moment  coefficient  slope  in  the  transonic  region  found  in  the  experiment  [5,  6], 
whereas  AP09  under-predicts  by  a  significant  margin. 

4.  Lift  Coefficient  at  Different  Speed  Regimes 

The  computed  lift  coefficients  at  different  speed  regimes  are  illustrated  in  Figure 
33,  34  and  35  for  AOA  of  0°,  2°  and  4°.  Both  ANSYS  CFX  and  AP09  correctly  predict 
the  upward  trend  of  Lift  Coefficient  with  increasing  AOA. 
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Figure  33.  Lift  Coefficient  at  Subsonic  Region. 


Figure  34.  Lift  Coefficient  at  Transonic  Region. 


43 


5.  Flow  Visualization 

Additional  insight  can  be  obtained  by  visualizing  the  flow  field  changes  with  the 
Mach  number  and  AOA.  Therefore,  Mach  contour  plots  in  the  XZ  plane  at  Mach 
Numbers  of  0.7,  0.9,  1  and  1.4  at  AOA  of  2°  using  IVM  are  presented  in  Figure  36,  37, 
38  and  39.  Figure  40,  41,  42  and  43  offer  an  alternate  perspective  in  the  Y-Z  plane.  From 
Figures  44  and  45,  one  can  observe  the  change  in  the  base  velocity  vectors  as  the  AOA  is 
changed  from  0°  and  2°. 


44 


Mach  Number 
Contour  1 


■ 

■ 


8029e-001 
7.607e-001 
7.184e-001 
6.761e-001 
6  339e-001 
5.916e-001 
5  494e-001 
5.071  e-001 
4.648e-001 
4  226e-001 
3.803e-001 
3  381e-001 
2.958e-001 
2.536e-001 
2  1 1 3e-001 
1.690e-001 
1  268e-001 
8  452e-002 
4.226e-002 
1  000e-015 


1.000 


Figure  36.  Mach  Contour  at  M=0.7,  AOA  of  2°,  IVM. 
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Figure  37.  Mach  Contour  at  M=0.9,  AOA  of  2°,  IVM. 
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Figure  38.  Mach  Contour  at  M=l,  AOA  of  2°,  IVM. 
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Mach  Contour  at  M=1 .4,  AOA  of  2°,  IVM. 
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Figure  39. 
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Figure  40. 
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Side  Profile-Mach  Contour  at  M=0.7, 
AOA  of  2°,  IVM. 
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Figure  41 .  Side  Profile-Mach  Contour  at  M=0.9, 

AOA  of  2°,  IVM. 
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Figure  42. 


Side  Profile-Mach  Contour  at  M=l, 
AOA  of  2°,  IVM. 
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Figure  43.  Side  Profile-Mach  Contour  at  M=1 .4, 

AOA  of  2°,  IVM. 
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Figure  44.  Base  Velocity  Vector  Flow  at  AOA  of  0°,  IVM. 


Figure  45.  Base  Velocity  Vector  Flow  at  AOA  of  2°,  IVM. 
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B. 


PART  II-COMPARISON  OF  AOA  GENERATION  METHODS 


In  Part  I,  the  simulations  were  perfonned  using  the  IVM  method  for  AOA 
Generation.  In  Part  II,  additional  simulations  for  the  Standard  Nose  Configuration  were 
perfonned  with  the  BIR  method.  Using  BIR,  the  body  of  interest  is  rotated  to  the  desired 
AOA.  Results  from  the  two  methods  are  compared  in  Figure  46,  47  and  48.  It  is  observed 
from  Figure  46  that  the  two  methods  produce  essentially  identical  total  drag  coefficients. 
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Figure  46.  Comparison  of  Total  Drag  Coefficient  for  IVM  and  BIR  AOA 

Generation  Methods. 


The  agreement  in  predicted  normal  coefficient  slope,  shown  in  Figure  47  and  in 
pitching  moment  coefficient  slope,  shown  in  Figure  48,  is  somewhat  less  satisfactory. 
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Figure  47.  Comparison  of  Normal  Force  Coefficient  Slope  for  IVM  and  BIR 

AOA  Methods. 
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Figure  48.  Comparison  of  Pitching  Moment  Coefficient  Slope  for  IVM  and  BIR 

AOA  Generation  Methods. 
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Mach  contours  at  Mach  Numbers  of  0.7,  0.9,  1  and  1.4  at  AOA  of  2°  using 
BIR  are  illustrated  in  Figures  49,  50,  51  and  52.  From  the  plots,  it  is  observed  that 
the  Mach  contours  of  BIR  are  in  good  agreement  with  that  of  IVM. 
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Figure  49.  Mach  Contour  at  M=0.7,  AOA  of  2°,  BIR. 
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Figure  50.  Mach  Contour  at  M=0.9,  AOA  of  2°,  BIR. 
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Figure  5 1 .  Mach  Contour  at  M=l,  AOA  of  2°,  BIR. 
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Figure  52. 


Mach  Contour  at  M=1.4,  AOA  of  2°,  BIR. 


1.  Analysis  of  Results 

IVM  is  believed  to  provide  a  better  prediction  than  BIR.  In  IVM,  the  free-stream 
flow  direction  is  deliberately  offset  from  the  X  axis  direction.  The  formula  stated  in 
Equation  2.3  and  2.4  is  applied  so  that  the  forces  calculated  are  adjusted  to  account  for 
the  AOA.  “force  y”  in  Function  Calculator  is  equal  to  nonnal  force  on  the  body  of 
interest.  In  BIR,  the  generated  “force  y”  is  the  lift  force  perpendicular  to  the  free-stream 
velocity.  Theoretically,  by  applying  the  formula  in  2.3  and  2.4,  the  result  should  be  the 
same.  However,  the  tilted  body  in  BIR  increased  the  meshing  complexity,  evident  in  the 
increase  of  mesh  nodes  and  elements  as  compared  to  IVM.  The  difference  in  mesh  profile 
is  believed  to  be  the  reason  for  the  minor  variation  between  the  two  methods. 

C.  PART  III-MODIFIED  NOSE  CONFIGURATION 

Modification  to  the  projectile  body  is  not  desirable  as  it  will  discontinue  the  usage 
of  existing  projectile  stockpile.  Thus,  modification  was  made  only  on  the  fuze.  It  is 
proposed  that  a  new  design  fuze  is  able  to  tilt  at  various  angles  to  allow  lift  generation. 
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Modification  to  the  fuze  model  was  perfonned  using  SolidWorks  as  illustrated  in 
Figure  53.  In  order  to  check  whether  lift  is  generated  by  tilting  the  nose  of  the  projectile, 
the  fuze  is  made  to  tilt  at  angles  of  2°,  4°,  6°,  8°  and  10°. 


Figure  53.  Modified  Fuze  Model. 


For  the  Modified  Nose  Configuration,  ANSYS  CFX  and  AP09  simulations  were 
perfonned  on  a  modified  M549  projectile  model  for  various  AOA  and  Mach  numbers. 
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1.  Lift  Coefficient  versus  Mach  Number  at  Nose  Tilt  Angles  of  0°,  2°,  4°, 
6°,  8°  and  10° 

The  computed  lift  coefficient  versus  Mach  number  distributions  are  illustrated  in 
Figure  54.  They  were  obtained  using  eight  data  points  from  Mach  0.7  to  Mach  1.4  with 
an  interval  of  0.1. 
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Figure  54.  Lift  Coefficient  versus  Mach  Number. 

From  Figure  54,  it  is  seen  that,  as  expected,  in  subsonic  and  transonic  flight,  lift 
increased  as  the  nose  tilt  was  increased  to  2°  and  4°,  but  then  started  to  diminish  with  a 
further  tilt  angle  increase.  In  contrast,  in  supersonic  flight,  the  tilted  nose  reduced  the 
lift.  An  explanation  of  the  results  shown  in  Figure  54  is  offered  in  Figure  55  and  56. 
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2.  Analysis  of  Results 

Velocity  vector  plots  at  nose  tilt  angles  of  0°,  4°  and  6°  at  Mach  0.7  are  shown  in 

Figure  55. 
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Figure  55.  Velocity  Vector  Plot  for  Nose  Tilt  Angles  of  0°,  4°  and  6°  at  Mach  0.7. 


It  is  seen  that  the  flow  is  attached  when  the  nose  tilt  angle  is  0°,  2°  and  4°. 
However  flow  separation  is  observed  at  nose  tilt  angle  of  6°  and  beyond.  Apparently  this 
flow  separation  prevented  further  improvement  in  lift  as  the  nose  tilt  angles  increased. 
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Figure  56.  Velocity  Vector  Plot  for  Nose  Tilt  Angles  of  0°,  4°  and  6°  at  Mach  1 .2. 


In  supersonic  flow  as  seen  from  Figure  56,  tilting  the  nose  caused  flow  separation 
for  all  tilt  angles  and  therefore  led  to  a  reduction  of  lift.  However,  further  analysis  is 
required  to  ascertain  the  numerical  accuracy  of  the  computations. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  key  objectives  of  this  thesis  were  successfully  completed.  The  non-spinning 
M549  projectile’s  stability  derivatives  were  successfully  validated  using  ANSYS  CFX 
and  AP09.  Comparison  studies  also  were  perfonned  using  two  different  AOA  generation 
methods  in  ANSYS  CFX.  Furthermore,  the  feasibility  of  lift  generation  via  projectile 
nose  modification  was  explored  using  ANSYS  CFX. 

In  general,  ANSYS  CFX  was  observed  to  provide  better  trend  prediction  than 
AP09.  In  this  thesis,  a  prediction  error  of  +/- 20%  was  observed.  This  was  due  to 
restricted  mesh  quality  imposed  by  computing  limits.  This  highlighted  the  need  to  have 
sufficient  computing  power  to  achieve  better  mesh  quality  for  accuracy.  The  IVM  AOA 
generation  method  provided  a  better  prediction  than  the  BIR  method.  Differences  in  mesh 
profiles  contributed  to  the  minor  variation  between  the  two  methods. 

Lift  generation  via  nose  tilt  modification  was  predicted  for  nose  tilt  angles  of  2°, 
4°,  6°,  8°  and  10°.  The  onset  of  flow  separation  was  found  to  impose  limits  on  this 
method  of  lift  generation. 

B.  RECOMMENDATIONS 

The  proposed  next  stage  of  investigation  is  to  simulate  a  spinning  M549  projectile. 
This  is  necessary  to  detennine  the  influence  of  the  spinning  motion  on  the  aerodynamics 
of  the  projectile. 

In  this  thesis,  it  was  observed  that  both  ANSYS  CFX  and  AP09  have  limitations 
in  providing  accurate  prediction.  Alternate  prediction  codes  should  be  considered  in  the 
next  stage  of  investigation.  Commercially  available  prediction  codes  like  Projectile 
Rocket  Ordnance  Design  and  Analysis  System  (PRODAS)  offer  more  than  the  standard 
aerodynamic  coefficients  prediction.  It  allows  users  to  conduct  launch  dynamic  and 
trajectories  simulation.  These  simulations  provide  essential  infonnation  for  the  study  of 
the  aerodynamic  properties  of  the  projectile. 
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APPENDIX  A.  MESH  INPUT  DATA 


The  following  mesh  input  data  is  specified  for  this  thesis.  Other  required  input 
data  not  mentioned  in  the  list  is  set  as  default. 


Sizing 

Use  Advanced 

On  Curvature 

Relevance  Center 

Fine 

Curvature  Nonnal  Angle 

4  Deg 

Min  Size 

le-5m 

Max  Face  Size 

0.1m 

Max  Size 

0.1m 

Growth  Rate 

1.2 

Min  Edge 

1.55e-4m 

Inflation 

Inflation  Option 

Total  Thickness 

Number  of  Layer 

40 

Growth  Rate 

1.1 

Max  Thickness 

le-3m 

Vertex  Sizing 

Element  Size 

le-2m 

Sphere  Radius 

0.4m 

Table  1.  Mesh  Input  Data. 


63 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


64 


APPENDIX  B.  ANSYS  CFX  COMMAND  LANGUAGE  FOR  RUN 

(LOW  MACH  NUMBER) 


LIBRARY: 

MATERIAL:  Air  Ideal  Gas 
Material  Description  =  Air  Ideal  Gas  (constant  Cp) 
Material  Group  =  Air  Data,  Calorically  Perfect  Ideal  Gases 
Option  =  Pure  Substance 
Thermodynamic  State  =  Gas 
PROPERTIES: 

Option  =  General  Material 
EQUATION  OF  STATE: 

Molar  Mass  =  28.96  [kg  kmolA-l] 

Option  =  Ideal  Gas 
SPECIFIC  HEAT  CAPACITY: 

Option  =  Value 

Specific  Heat  Capacity  =  1.0044E+03  [J  kgA-l  KA-1] 
Specific  Heat  Type  =  Constant  Pressure 
REFERENCE  STATE: 

Option  =  Specified  Point 
Reference  Pressure  =  1  [atm] 

Reference  Specific  Enthalpy  =  0.  [J/kg] 

Reference  Specific  Entropy  =  0.  [J/kg/K] 

Reference  Temperature  =  25  [C] 

DYNAMIC  VISCOSITY: 

Dynamic  Viscosity  =  1.831E-05  [kgmA-l  sA-l] 

Option  =  Value 

THERMAL  CONDUCTIVITY: 

Option  =  Value 

Thermal  Conductivity  =  2.61E-2  [W  mA-l  KA-1] 
ABSORPTION  COEFFICIENT: 

Absorption  Coefficient  =  0.01  [mA-l] 

Option  =  Value 

SCATTERING  COEFFICIENT: 

Option  =  Value 

Scattering  Coefficient  =  0.0  [mA-l] 

REFRACTIVE  INDEX: 

Option  =  Value 

Refractive  Index  =  1.0  [m  mA-l] 

FLOW:  Flow  Analysis  1 
SOLUTION  UNITS: 
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Angle  Units  =  [rad] 

Length  Units  =  [m] 

Mass  Units  =  [kg] 

Solid  Angle  Units  =  [sr] 

Temperature  Units  =  [K] 

Time  Units  =  [s] 

ANALYSIS  TYPE: 

Option  =  Steady  State 
EXTERNAL  SOLVER  COUPLING: 

Option  =  None 
DOMAIN :  Default  Domain 
Coord  Frame  =  Coord  0 
Domain  Type  =  Fluid 
Location  =  B36 

BOUNDARY :  Default  Domain  Default 
Boundary  Type  =  WALL 
Location  =  F87.36,F88.36,F89.36,F90.36,F91.36 
BOUNDARY  CONDITIONS: 

HEAT  TRANSFER: 

Option  =  Adiabatic 
MASS  AND  MOMENTUM: 

Option  =  No  Slip  Wall 
WALL  ROUGHNESS: 

Option  =  Smooth  Wall 
BOUNDARY :  bottom 
Boundary  Type  =  INLET 
Location  =  bottom 
BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
HEAT  TRANSFER: 

Option  =  Static  Temperature 
Static  Temperature  =  288.15  [K] 

MASS  AND  MOMENTUM: 

Option  =  Cartesian  Velocity  Components 
U  =  -238  [m  sA-l] 

V  =  8.33  [m  sA-l] 

W  =  0  [m  sA-l] 

TURBULENCE: 

Option  =  Medium  Intensity  and  Eddy  Viscosity  Ratio 
BOUNDARY:  inlet 
Boundary  Type  =  INLET 
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Location  =  inlet 
BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
HEAT  TRANSFER: 

Option  =  Static  Temperature 
Static  Temperature  =  288.15  [K] 

MASS  AND  MOMENTUM: 

Option  =  Cartesian  Velocity  Components 
U  =  -238  [m  sA-l] 

V  =  8.33  [m  sA-l] 

W  =  0  [m  sA- 1  ] 

TURBULENCE: 

Option  =  Medium  Intensity  and  Eddy  Viscosity  Ratio 
BOUNDARY:  outlet 
Boundary  Type  =  OUTLET 
Location  =  outlet 
BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
MASS  AND  MOMENTUM: 

Option  =  Static  Pressure 
Relative  Pressure  =  0  [Pa] 

BOUNDARY:  syml 
Boundary  Type  =  SYMMETRY 
Location  =  syml 
BOUNDARY:  sym2 
Boundary  Type  =  WALL 
Location  =  sym2 
BOUNDARY  CONDITIONS: 

HEAT  TRANSFER: 

Option  =  Adiabatic 
MASS  AND  MOMENTUM: 

Option  =  No  Slip  Wall 
WALL  ROUGHNESS: 

Option  =  Smooth  Wall 
BOUNDARY:  top 
Boundary  Type  =  OUTLET 
Location  =  top 

BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
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MASS  AND  MOMENTUM: 

Option  =  Average  Static  Pressure 
Pressure  Profile  Blend  =  0.05 
Relative  Pressure  =  0  [Pa] 

PRESSURE  AVERAGING: 

Option  =  Average  Over  Whole  Outlet 
DOMAIN  MODELS: 

BUOYANCY  MODEL: 

Option  =  Non  Buoyant 
DOMAIN  MOTION: 

Option  =  Stationary 
MESH  DEFORMATION: 

Option  =  None 
REFERENCE  PRESSURE: 

Reference  Pressure  =  1  [atm] 

FLUID  DEFINITION:  Fluid  1 
Material  =  Air  Ideal  Gas 
Option  =  Material  Library 
MORPHOLOGY: 

Option  =  Continuous  Fluid 
FLUID  MODELS: 

COMBUSTION  MODEL: 

Option  =  None 

HEAT  TRANSFER  MODEL: 

Include  Viscous  Work  Term  =  On 
Option  =  Total  Energy 
THERMAL  RADIATION  MODEL: 

Option  =  None 
TURBULENCE  MODEL: 

Option  =  SST 

TRANSITIONAL  TURBULENCE: 

Option  =  Gamma  Theta  Model 
TRANSITION  ONSET  CORRELATION: 
Option  =  Langtry  Menter 
TURBULENT  WALL  FUNCTIONS: 

High  Speed  Model  =  Off 
Option  =  Automatic 
INITIALISATION: 

Option  =  Automatic 
INITIAL  CONDITIONS: 

Velocity  Type  =  Cartesian 
CARTESIAN  VELOCITY  COMPONENTS: 
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Option  =  Automatic  with  Value 
U  =  -238  [m  sA-l] 

V  =  8.33  [m  sA-l] 

W  =  0  [m  sA-l] 

STATIC  PRESSURE: 

Option  =  Automatic 
TEMPERATURE: 

Option  =  Automatic 

TURBULENCE  INITIAL  CONDITIONS: 

Option  =  Medium  Intensity  and  Eddy  Viscosity  Ratio 
OUTPUT  CONTROL: 

RESULTS: 

File  Compression  Level  =  Default 
Option  =  Standard 
SOLVER  CONTROL: 

Turbulence  Numerics  =  High  Resolution 
ADVECTION  SCHEME: 

Option  =  High  Resolution 
COMPRESSIBILITY  CONTROL: 

High  Speed  Numerics  =  On 
CONVERGENCE  CONTROL: 

Length  Scale  Option  =  Conservative 
Maximum  Number  of  Iterations  =100 
Minimum  Number  of  Iterations  =  1 
Timescale  Control  =  Auto  Timescale 
Timescale  Factor  =  10 
CONVERGENCE  CRITERIA: 

Residual  Target  =  le-07 
Residual  Type  =  RMS 
DYNAMIC  MODEL  CONTROL: 

Global  Dynamic  Model  Control  =  On 
EXPERT  PARAMETERS: 
max  continuity  loops  =  3 
topology  estimate  factor  =1.1 
COMMAND  FIT  E: 

Version  =  13.0 
Results  Version  =  13.0 
SIMULATION  CONTROL: 

EXECUTION  CONTROL: 

EXECUTABLE  SELECTION: 

Double  Precision  =  On 
INTERPOLATOR  STEP  CONTROL: 
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Runtime  Priority  =  Standard 
PARALLEL  HOST  LIBRARY: 

HOST  DEFINITION:  skypc 
Remote  Host  Name  =  SKY-PC 
Host  Architecture  String  =  winnt-amd64 
Installation  Root  =  C:\Program  Files\ANSYS  Inc\v%v\CFX 
PARTITIONER  STEP  CONTROL: 

Multidomain  Option  =  Independent  Partitioning 
Runtime  Priority  =  Standard 
EXECUTABLE  SELECTION: 

Use  Large  Problem  Partitioner  =  Off 
PARTITIONING  TYPE: 

MeTiS  Type  =  k-way 
Option  =  MeTiS 
Partition  Size  Rule  =  Automatic 
Partition  Weight  Factors  =  0.50000,  0.50000 
RUN  DEFINITION: 

Run  Mode  =  Full 

Solver  Input  File  =  Fluid  Flow  CFX.def 
SOLVER  STEP  CONTROL: 

Runtime  Priority  =  Standard 
PARALLEL  ENVIRONMENT: 

Number  of  Processes  =  2 

Start  Method  =  HP  MPI  Local  Parallel 

Parallel  Host  List  =  skypc*2 
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APPENDIX  C.  ANSYS  CFX  COMMAND  LANGUAGE  FOR  RUN 

(HIGH  MACH  NUMBER) 


LIBRARY: 

MATERIAL:  Air  Ideal  Gas 
Material  Description  =  Air  Ideal  Gas  (constant  Cp) 
Material  Group  =  Air  Data,  Calorically  Perfect  Ideal  Gases 
Option  =  Pure  Substance 
Thermodynamic  State  =  Gas 
PROPERTIES: 

Option  =  General  Material 
EQUATION  OF  STATE: 

Molar  Mass  =  28.96  [kg  kmolA-l] 

Option  =  Ideal  Gas 
SPECIFIC  HEAT  CAPACITY: 

Option  =  Value 

Specific  Heat  Capacity  =  1.0044E+03  [J  kgA-l  KA-1] 
Specific  Heat  Type  =  Constant  Pressure 
REFERENCE  STATE: 

Option  =  Specified  Point 
Reference  Pressure  =  1  [atm] 

Reference  Specific  Enthalpy  =  0.  [J/kg] 

Reference  Specific  Entropy  =  0.  [J/kg/K] 

Reference  Temperature  =  25  [C] 

DYNAMIC  VISCOSITY: 

Dynamic  Viscosity  =  1.831E-05  [kgmA-l  sA-l] 

Option  =  Value 

THERMAL  CONDUCTIVITY: 

Option  =  Value 

Thermal  Conductivity  =  2.61E-2  [W  mA-l  KA-1] 
ABSORPTION  COEFFICIENT: 

Absorption  Coefficient  =  0.01  [mA-l] 

Option  =  Value 

SCATTERING  COEFFICIENT: 

Option  =  Value 

Scattering  Coefficient  =  0.0  [mA-l] 

REFRACTIVE  INDEX: 

Option  =  Value 

Refractive  Index  =  1.0  [m  mA-l] 

FLOW:  Flow  Analysis  1 
SOLUTION  UNITS: 
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Angle  Units  =  [rad] 

Length  Units  =  [m] 

Mass  Units  =  [kg] 

Solid  Angle  Units  =  [sr] 

Temperature  Units  =  [K] 

Time  Units  =  [s] 

ANALYSIS  TYPE: 

Option  =  Steady  State 
EXTERNAL  SOLVER  COUPLING: 

Option  =  None 
DOMAIN :  Default  Domain 
Coord  Frame  =  Coord  0 
Domain  Type  =  Fluid 
Location  =  B36 

BOUNDARY :  Default  Domain  Default 
Boundary  Type  =  WALL 
Location  =  F87.36,F88.36,F89.36,F90.36,F91.36 
BOUNDARY  CONDITIONS: 

HEAT  TRANSFER: 

Option  =  Adiabatic 
MASS  AND  MOMENTUM: 

Option  =  No  Slip  Wall 
WALL  ROUGHNESS: 

Option  =  Smooth  Wall 
BOUNDARY :  bottom 
Boundary  Type  =  INLET 
Location  =  bottom 
BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
HEAT  TRANSFER: 

Option  =  Static  Temperature 
Static  Temperature  =  288.15  [K] 

MASS  AND  MOMENTUM: 

Option  =  Cartesian  Velocity  Components 
U  =  -476.13  [m  sA-l] 

V  =  16.67  [m  sA-l] 

W  =  0  [m  sA-l] 

TURBULENCE: 

Option  =  Medium  Intensity  and  Eddy  Viscosity  Ratio 
BOUNDARY:  inlet 
Boundary  Type  =  INLET 
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Location  =  inlet 


BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
HEAT  TRANSFER: 

Option  =  Static  Temperature 
Static  Temperature  =  288.15  [K] 

MASS  AND  MOMENTUM: 

Option  =  Cartesian  Velocity  Components 
U  =  -476.13  [m  sA-l] 

V  =  16.67  [m  sA-l] 

W  =  0  [m  sA- 1  ] 

TURBULENCE: 

Option  =  Medium  Intensity  and  Eddy  Viscosity  Ratio 
BOUNDARY:  outlet 
Boundary  Type  =  OUTLET 
Location  =  outlet 
BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
MASS  AND  MOMENTUM: 

Option  =  Static  Pressure 
Relative  Pressure  =  0  [Pa] 

BOUNDARY:  syml 
Boundary  Type  =  SYMMETRY 
Location  =  syml 
BOUNDARY:  sym2 
Boundary  Type  =  WALL 
Location  =  sym2 
BOUNDARY  CONDITIONS: 

HEAT  TRANSFER: 

Option  =  Adiabatic 
MASS  AND  MOMENTUM: 

Option  =  No  Slip  Wall 
WALL  ROUGHNESS: 

Option  =  Smooth  Wall 
BOUNDARY:  top 
Boundary  Type  =  OUTLET 
Location  =  top 

BOUNDARY  CONDITIONS: 

FLOW  REGIME: 

Option  =  Subsonic 
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MASS  AND  MOMENTUM: 

Option  =  Average  Static  Pressure 
Pressure  Profile  Blend  =  0.05 
Relative  Pressure  =  0  [Pa] 

PRESSURE  AVERAGING: 

Option  =  Average  Over  Whole  Outlet 
DOMAIN  MODELS: 

BUOYANCY  MODEL: 

Option  =  Non  Buoyant 
DOMAIN  MOTION: 

Option  =  Stationary 
MESH  DEFORMATION: 

Option  =  None 
REFERENCE  PRESSURE: 

Reference  Pressure  =  1  [atm] 

FLUID  DEFINITION:  Fluid  1 
Material  =  Air  Ideal  Gas 
Option  =  Material  Library 
MORPHOLOGY: 

Option  =  Continuous  Fluid 
FLUID  MODELS: 

COMBUSTION  MODEL: 

Option  =  None 

HEAT  TRANSFER  MODEL: 

Include  Viscous  Work  Term  =  On 
Option  =  Total  Energy 
THERMAL  RADIATION  MODEL: 

Option  =  None 
TURBULENCE  MODEL: 

Option  =  SST 

TRANSITIONAL  TURBULENCE: 

Option  =  Gamma  Theta  Model 
TRANSITION  ONSET  CORRELATION: 
Option  =  Langtry  Menter 
TURBULENT  WALL  FUNCTIONS: 

High  Speed  Model  =  Off 
Option  =  Automatic 
INITIALISATION: 

Option  =  Automatic 
INITIAL  CONDITIONS: 

Velocity  Type  =  Cartesian 
CARTESIAN  VELOCITY  COMPONENTS: 
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Option  =  Automatic  with  Value 
U  =  -476.13  [m  sA-l] 

V  =  16.67  [m  sA-l] 

W  =  0  [m  sA-l] 

STATIC  PRESSURE: 

Option  =  Automatic 
TEMPERATURE: 

Option  =  Automatic 

TURBULENCE  INITIAL  CONDITIONS: 

Option  =  Medium  Intensity  and  Eddy  Viscosity  Ratio 
OUTPUT  CONTROL: 

RESULTS: 

File  Compression  Level  =  Default 
Option  =  Standard 
SOLVER  CONTROL: 

Turbulence  Numerics  =  High  Resolution 
ADVECTION  SCHEME: 

Option  =  High  Resolution 
COMPRESSIBILITY  CONTROL: 

High  Speed  Numerics  =  On 
CONVERGENCE  CONTROL: 

Length  Scale  Option  =  Conservative 
Maximum  Number  of  Iterations  =  20 
Minimum  Number  of  Iterations  =  1 
Timescale  Control  =  Auto  Timescale 
Timescale  Factor  =  10 
CONVERGENCE  CRITERIA: 

Residual  Target  =  le-07 
Residual  Type  =  RMS 
DYNAMIC  MODEL  CONTROL: 

Global  Dynamic  Model  Control  =  On 
EXPERT  PARAMETERS: 
max  continuity  loops  =  3 
topology  estimate  factor  =1.1 
COMMAND  FIT  E: 

Version  =  13.0 
Results  Version  =  13.0 
SIMULATION  CONTROL: 

EXECUTION  CONTROL: 

EXECUTABLE  SELECTION: 

Double  Precision  =  On 
INTERPOLATOR  STEP  CONTROL: 
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Runtime  Priority  =  Standard 
PARALLEL  HOST  LIBRARY: 

HOST  DEFINITION:  skypc 
Remote  Host  Name  =  SKY-PC 
Host  Architecture  String  =  winnt-amd64 
Installation  Root  =  C:\Program  Files\ANSYS  Inc\v%v\CFX 
PARTITIONER  STEP  CONTROL: 

Multidomain  Option  =  Independent  Partitioning 
Runtime  Priority  =  Standard 
EXECUTABLE  SELECTION: 

Use  Large  Problem  Partitioner  =  Off 
PARTITIONING  TYPE: 

MeTiS  Type  =  k-way 
Option  =  MeTiS 
Partition  Size  Rule  =  Automatic 
Partition  Weight  Factors  =  0.50000,  0.50000 
RUN  DEFINITION: 

Run  Mode  =  Full 

Solver  Input  File  =  Fluid  Flow  CFX.def 
SOLVER  STEP  CONTROL: 

Runtime  Priority  =  Standard 
PARALLEL  ENVIRONMENT: 

Number  of  Processes  =  2 

Start  Method  =  HP  MPI  Local  Parallel 

Parallel  Host  List  =  skypc*2 
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APPENDIX  D.  OUTPUT  PLOTS  FOR  AOA  OF  0° 


Mach  Number 
Contour  1 


8.103e-001 
7  677e-001 

7  250e-001 
6.824e-001 
6  397e-001 
5.971  e-001 
5  544e-001 
5.118e-001 
4.691  e-001 
4  265e-001 
3.838e-001 
3.412e-001 
2.985e-001 
2  559e-001 
2  132e-001 
1.706e-001 
1.279e-001 

8  529e-002 
4.265e-002 
1  000e-015 


2.000  (m) 


Figure  57.  Mach  Contour  at  M=0.7,  AOA  of  0°,  IVM. 


Mach  Number 
Contour  1 


n 


9.602e-001 
9  096e-001 
8.591  e-001 
8.086e-001 
7  580e-001 
7.075e-001 
6  570e-001 
6.064e-001 
5.559e-001 
5  054e-001 
4.548e-001 
4.043e-001 
3.537e-001 
3.032e-001 
2  527e-001 
2.021e-001 
1.516e-001 
1.011  e-001 
5.054e-002 
1  000e-015 


C37 

♦  If 


Y 


0.500  1.500 


Figure  58.  Mach  Contour  at  M=0.8,  AOA  of  0°,  IVM. 
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Mach  Number 
Contour  1 


n 


1  138e+000 
1  079e+000 

1  019e+000 
9.587e-001 
8  988e-001 
8.388e-001 
7  789e-001 
7.190e-001 
6  591e-001 
5  992e-001 
5.393e-001 
4  793e-001 
4  194e-001 
3.595e-001 

2  996e-001 
2.397e-001 
1  798e-001 
1  198e-001 
5.992e-002 
1  000e-015 


if 


0.500  1.500 


Figure  59.  Mach  Contour  at  M=0.9,  AOA  of  0°,  IVM. 


Mach  Number 
Contour  1 


1  206e+000 
1  143e+000 
1  079e+000 
1.016e+000 
9  523e-001 
8  888e-001 
8.253e-001 
7.618e-001 
6.984e-001 
6.349e-001 
5.714e-001 
5.079e-001 
4  444e-00l 
3.809e-001 
3  174e-001 
2.539e-001 
1.905e-001 
1.270e-001 
6.349e-002 
1  000e-015 


2.000  (m) 


Figure  60. 


Mach  Contour  at  M=l,  AOA  of  0°,  IVM. 
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n 


Mach  Number 

Contour  1 


1  256e+000 
1.190e+000 
1.124e+000 
1.058e+000 
9  915e-001 
9.254e-001 
8  593e-001 
7  932e-001 
7.271e-001 
6  610e-001 
5  949e-001 
5  288e-001 
4.627e-001 
3.966e-001 
3  305e-001 

2  644e-001 
1.9836-001 
1  322e-001 
6.610e-002 
1  OOOe-015 


Figure  61.  Mach  Contour  at  M=l.l,  AOA  of  0°,  IVM. 


Mach  Number 
Contour  1 


1.346e+000 

1  275e+000 
1.205e+000 
1 . 1 34e+000 
1.063e+000 
9.920e-001 
9.212e-001 
8  503e-001 
7.794e-001 
7  086e-001 
6.377e-001 
5.669e-001 
4.960e-001 
4.252e-001 
3.543e-001 
2.834e-001 

2  126e-001 
1.417e-001 
7.086e-002 
1, 000e-015 


Figure  62. 


Mach  Contour  at  M=1.2,  AOA  of  0°,  IVM. 
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APPENDIX  E.  OUTPUT  PLOTS  FOR  AOA  OF  2° 


fl 


Mach  Number 
Contour  1 


8.029e-001 
7  607e-001 
7  184e-001 
6.761e-001 
6  339e-001 
5.916e-001 
5.494e-001 
5.071  e-001 
4.648e-001 
4  226e-001 
3.803e-001 
3.381e-001 
2.958e-001 
2.536e-001 
2  1 1 3e-001 
1  690e-001 
1  268e-001 
8.452e-002 
4  226e-002 
1  000e-015 


Figure  63.  Mach  Contour  at  M=0.7,  AOA  of  2°,  IVM. 


Mach  Number 
Contour  1 


9  546e-001 
9  043e-001 
8.541e-001 
8.039e-001 
7  536e-001 
7  034e-001 
6  531e-001 
6.029e-001 
5  527e-001 
5  024e-001 
4.522e-001 
4  019e-001 
3.517e-001 
3.014e-001 
2  512e-001 
2.010e-001 
1  507e-001 
1  005e-001 
5.024e-002 
1  000e-015 


2.000  (m) 


0.500 


Figure  64.  Mach  Contour  at  M=0.8,  AOA  of  2°,  IVM. 
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Mach  Number 
Contour  1 


1  138e+000 
1.078e+000 
1.019e+000 
9.587e-001 
8  987e-001 
8.388e-001 
7.789e-001 
7.190e-001 
6.591  e-001 
5  992e-001 
5.392e-001 
4  793e-001 
4  194e-001 
3.595e-001 

2  996e-001 
2.397e-001 
1  797e-001 
1.198e-001 
5.992e-002 
1  000e-015 


Figure  65. 


1.000 


2.000  (m) 


0.500  1.500 

Mach  Contour  at  M=0.9,  AOA  of  2°,  IVM. 


Mach  Number 
Contour  1 


Ft 


1  210e+000 
1  146e+000 
1  083e+000 
1.019e+000 
9  552e-001 
8.915e-001 
8.278e-001 
7.642e-001 
7.005e-001 
6  368e-001 
5.731  e-001 
5.094e-001 
4  458e-001 
3.821  e-001 
3  184e-001 
2.547e-001 
1  910e-001 
1  274e-001 
6.368e-002 
1  000e-015 


Figure  66 


Mach  Contour  at  M=l,  AOA  of  2°,  IVM. 
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Mach  Number 
Contour  1 


1.278e+00Q 
1  211e+000 
1.143e+000 
1.076e+000 
1  009e+000 
9  416e-001 
8  743e-001 
8  070e-001 
7  398e-001 
6  725e-001 
6  053e-001 
5  380e-001 
4  708e-001 
4.035e-001 
3  363e*001 
2.690e-001 
2.018e-001 
1  345e-001 
6.725e-002 
1  OOOe-015 


Figure  67.  Mach  Contour  at  M=l.l,  AOA  of  2°,  IVM. 


Mach  Number 
Contour  1 


1.353e+000 
1  282e+000 
1.210e+000 
1.139e+000 

1  068e+000 
9.968e-001 
9  256e-001 
8.544e-001 
7.832e-001 
7  120e-001 
6.408e-001 
5  696e-001 
4  984e-00l 
4.272e-001 
3  560e-001 
2.848e-001 

2  136e-001 
1.424e-001 
7.120e-002 
1  000e-015 


Figure  68. 


Mach  Contour  at  M=1.2,  AOA  of  2°,  IVM. 
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Mach  Number 
Contour  1 


1.445e+000 
1  369e+000 
1  293e+000 
1  217e+000 
1  141e+000 
1  064e+000 
9  885e-001 
9  124e-001 
8.364e-001 
7  604e-001 
6  843e-001 

6  083e-001 
5  322e-001 
4  562e-001 
3  802e-001 
3.041e-001 
2.281e-001 
1  521e-001 

7  604e-002 
1  000e-015 


Figure  69. 


Mach  Contour  at  M=1 .3  AOA  of  2°,  IVM. 


Mach  Number 
Contour  1 


1.538e+000 
1  457e+000 
1  376e+000 
1.295e-K)00 
1  215e+000 
1 . 1 34e*000 
1.053e+000 
9.716e-001 
8.906e-001 
8.097e-001 
7.287e-001 
6  477e-001 
5.668e-00t 
4.858e-001 
4  048e-001 
3.239e-001 
2.429e-001 
1  619e-001 
8.097e-002 
1  000e-015 


Figure  70. 


Mach  Contour  at  M=1 .4  AOA  of  2°,  IVM. 


84 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Professor  and  Chainnan  Knox  T.  Millsaps 

Department  of  Mechanical  Engineering  and  Aeronautical  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

4.  Professor  Maximilian  F.  Platzer 

Department  of  Mechanical  Engineering  and  Aeronautical  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

5.  Professor  Anthony  J.  Gannon 

Department  of  Mechanical  Engineering  and  Aeronautical  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

6.  Professor  Yeo  Tat  Soon 

Temasek  Defence  Systems  Institute 
National  University  of  Singapore 
tdshead@nus  .edu.  sg 

7.  Ms.  Tan  Lai  Poh 

Temasek  Defence  Systems  Institute 
National  University  of  Singapore 


85 


